捷联惯性导航中一种算法的推证

捷联惯性导航系统中,常常用到解算四元数矩阵微分方程的问题。为此给出了解四元数矩阵微分方程的一个强有力的算式，并给出了详细的数学推导和证明。     

基于最优控制的四元数据误差传播捷联矩阵算法分析

捷联系统中的比力、角速率的量测值都是沿船舶坐标系的，如欲求得船舶位置和姿态就应知道其在当地地理坐标系中的投影。无论从计算速度还是从计算精度上考虑，用四元数法得到捷联短阵比其他方法更为有利。该就四元数计算所引起的误差及其传播方式与四元数向量的正交计算进行了讨论，从而获得了一种捷联矩阵算法。     

弹塑性板壳结构非线性有限元分析

本文根据塑性流动理论的基本公式，由隐式积分导出了与路径无关的变量更新算法和一致切线模量。采用单元广义应力应变直接离散塑性流动定律，构造了杂交应力单元一致切线刚度矩阵的显式表达式，编制了结构有限元程序ＳＡＦＥ，数值算例表明：本文的计算方法和计算程序是正确可靠的，可用于弹塑性板壳结构的非线性分析，计算结果屈曲临界载荷和极限承载能力。     

基于区域分解算法的地下水耦合模型及其应用

根据岩体中不同区域裂隙发育规模的差异分为连续区域和离散区域,在不同区域中使用不同的地下水运动数学模型,应用区域分解算法来解决这类问题。其中连续区域采用了等效连续介质模型,离散区域采用了随机裂隙网络模型,通过区域公共边界上水位和流量连续的条件将两模型耦合求解。将基于区域分解算法的耦合模型应用于锦屏水电站坝址区三维渗流场的模拟中,通过钻孔观测水位和计算水位的对比发现,该方法是有效的,能够应用于实际工程。     

关于激光捷联惯性导航系统姿态算法的分析

圆锥误差和量化误差是激光捷联惯性导航系统姿态解算误差的两个最主要的误差源.从分析圆锥误差产生的机理出发,分别分析了以角度和角速度为计算参数的圆锥误差补偿算法,并对量化误差对圆锥误差补偿算法的影响进行了研究.通过理论分析和数字仿真,得出在实际工程应用中,采用角速度为输入信息的激光捷联惯性导航系统姿态算法应该在考虑量化误差的情况下,采用以角速度为计算参数的圆锥误差补偿算法.     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

局部人工边界稳定性的一种分析方法

采用时域局部人工边界同有限元或有限差分相结合的方法分析无界域中的波动问题时可能出现由边界引起的数值失稳现象．本文提出一种分析局部人工边界稳定性的实用方法．该方法从时域逐步积分格式出发，考虑了边界节点运动和内部节点运动的耦合以及离散化的影响，可以给出局部人工边界在多维波动数值模拟中的稳定性近似准则．文中对多次透射边界在一维及二维出平面波动模拟中稳定性的分析结果，表明了该方法的合理性     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

延性单晶非均匀变形的三维有限元模拟

对延性单晶在拉伸载荷作用下的应变局域化和颈缩等非均匀变形过程进行了三维有限元数值模拟。将相关晶体塑性本构模型及一种新的数值积分方法补充到ABAQUS6.1商用有限元软件中。该方法的特点是，利用晶体塑性的动力学方程，获得一个关于晶体弹性变形梯度的演化方程，采用半隐式积分方案进行求解。本文推导出一种新的应力变本构矩阵。按此方式更新本构矩阵，计算速度和计算稳定性大大提高。加载方式，边界条件和变形程度等因素影响着滑移系的启动状况，这是平面模型所不能预测的。本文利用三维有限元方法模拟了不同取向下滑移系的启动状况，全面地考虑了FCC单晶材料12个可能滑移系在变形过程中的启动状况，合理地模拟了FCC面心立方单晶沿不同取向加载时晶轴旋转导致的应变局域化和颈缩等非均匀变形过程。     

A Non-linear Algorithm of Photoelastic Tomography for the Axisymmetric Problem

A non-linear algorithm of photoelastic tomography for the measurement of axisymmetric stress fields has been elaborated. It is free of any assumptions concerning the value of the birefringence or rotation of the principal stress axes along the light rays. The algorithm is based on the measurement of characteristic directions and phase retardation in two parallel sections of the test object. Stress components are presented in the form of power series along the radial coordinate. A differential evolution algorithm has been used for finding the stress field parameters, which fit the measurement data best. Application of the method is illustrated by residual stress measurement in a drinking glass.